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INTERACTION  OF  ELECTROMAGNETIC  WAVES 
WITH  A MOVING  IONIZATION  FRONT 

I.  Introduction 

The  refractive  index  non-uniformity  at  the  edge  of  a 
plasma  can  act  as  a mirror,  reflecting  incident  electromagnetic  waves. 
If  the  plasma  edge  is  moving  at  velocity  U,  the  reflected  wave  will  be 
doubly  Doppler  shifted,1  so  that  its  frequency2  oj*  is  related  to  the 
incident  frequency  x*  by 


(D*  1 + 3 cos  0*  . 

-S-- i-,  (1) 

u#  1 - B cos  9* 

where  S = U/c  and  0*,  9*  are  the  angles  of  incidence  and  reflectioi  3 
Furthermore,  a similar  Doppler  relation  occurs  between  the  duration  of 
the  reflected  pulse  t*  and  that  of  the  incident  pulse  t*, 


(2) 


Equation  (2)  applies  in  free  space;  the  situation  in  a waveguide  is 
somewhat  different  and  will  be  discussed  later.  Thus  an  upshifted 
wave  pulse  will  also  be  compressed,  which  tends  to  increase  the  peak 
reflected  power. 

Granatstein,  et  al.,3  in  a recent  experiment,  have  taken  advantage 
of  this  phenomenon  to  produce  intense  pulses  of  8mm  radiation  by 
reflecting  an  incident  3.2  cm  microwave  pulse  off  the  front  of  a 
high-current  relativistic  electron  beam,  which  plays  the  role  of  the 
overdense  plasma.4  This  experiment  is  part  of  a recent  surge  of 


Note.  Manuscript  submitted  May  23.  1977. 
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interest  in  the  exploitation  of  advances  in  relativistic  electron 
beam  technology  to  generate  intense  millimeter  and  submillimeter 
radiation.5  However,  we  pointed  out  recently6  that  the  plasma  itself 
need  not  be  moving,  for  the  wave  reflected  off  the  plasma  edge  to  be 
Doppler  shifted.  It  is  only  necessary  that  the  boundary  between  the 
plasma  and  the  surrounding  neutral  medium  be  moving.  Thus  an  alternative 
method  for  upshifting  a microwave  beam  is  to  reflect  it  from  a moving 
ionization  or  recombination  front  in  a stationary  gas,  i.e.  a moving 
interface  between  regions  where  the  gas  ( stationary  in  the  laboratory 
frame)  is  unionized  and  ionized.  (In  the  case  of  an  ionization  front, 
the  interface  moves  toward  the  neutral  gas,  U > 0;  a recombination 
front  moves  toward  the  plasma,  U < 0.)  In  Ref.  6 we  proposed  that 
this  might  be  a practical  method  of  generating  millimeter  and  sub- 
millimeter radiation  and  presented  preliminary  theoretical  calculations 
of  the  process.  In  the  present  paper,  these  are  extended  to  a detailed 
and  comprehensive  theory.  Throughout  this  paper,  the  plasma  is  treated 
as  a cold  and  collisionless  fluid.  Also  the  theory  is  linear  in  the 
amplitude  of  the  electromagnetic  wave,  so  that  no  scattering  processes 
other  than  refractive  index  variation  are  considered.  Except  for  these 
assumptions,  the  theory  is  quite  general,  and  no  mathematical  approxi- 
mations are  needed. 

In  Section  II  we  treat  the  case  of  a wave  incident  normally  on  a 
sharp,  infinite  planar,  moving  ionization  or  recombination  front  (i.e. 
a discontinuity  in  the  electron  density).  This  is  the  simplest  case 
to  understand.  In  Section  IIA  the  reflection  coefficient  is  calculated 
exactly.  It  is  demonstrated  that  the  Doppler  shift  and  reflected  pulse 


2 


duration  are  given  by  Eqs.  (l)  and  (2)  with  S*  = 9*  = 0,  the  same  as 
for  a moving  plasma,3  but  in  other  respects  reflection  off  a moving 
ionization  front  differs  significantly  from  reflection  off  a moving 
plasma.  A moving  overdense  plasma  (the  underdense  case  is  also  treated 
in  Sec.  II,  but  will  not  be  discussed  in  this  introduction),  acts  as  a 
perfect  moving  mirror,  reflecting  and  upshifting  each  photon.  Thus, 
according  to  Eqs.  (l)  and  (2),  the  energy  e*  and  power  of  the 
reflected  pulse  are  related  to  ihe  incident  e*  and  by 


£*/£*  = 

(3) 

P^/P*  = (t^/uS*)2 

(M 

Equations  (3)  or  (4)  hold  for  either  an  oncoming  or  receding  "mirror". 

For  an  oncoming  mirror  e*  exceeds  £*,  with  the  extra  energy  being  supplied  by  the 

kinetic  energy  of  the  moving  plasma;  for  the  receding  mirror,  energy  is 

absorbed  by  the  recoiling  mirror.  However,  a moving  ionization  front 

in  a stationary  plasma  has  no  kinetic  energy  to  share  with  the  wave; 

consequently,  one  must  have  £*  ^ e*.  The  complete  calculation  does 

indeed  show  that  this  is  true.  Furthermore,  the  physical  processes 

at  work  are  shown  to  be  different  for  the  case  of  an  oncoming  ionization 

front  (U  > 0,  0 > 0),  as  opposed  to  a receding  recombination  front 

(U  <0,  S < 0).  Therefore  £*/£*  and  P^/P*  cannot  be  represented  by  a 

single  analytic  expression  valid  for  both  U > 0 and  U < 0,  as  in  Eqs. 

1-4).  We  find  that  Eqs.  (3,4)  hold  for  a recombination  front  but 
that  for  an  ionization  front. 
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€*/e*  = a* /cu*. 

(5) 

P*/P*  * 1. 
r i 

(6) 

The  difference  between  the  cases  U < 0 and  0 > 0 is  explained  in  terms 
of  the  excitation  in  the  latter  case  of  a third  normal  mode  in  the 
plasma,  in  addition  to  the  usual  evanescent  and  growing  electromagnetic 
modes.  The  existence  of  this  third  wave  has  not  been  previously 
recognized . 

In  Section  IIB  it  is  demonstrated  explicitly  that  the  reflection 
and  transmission  coefficients  satisfy  energy  conservation.  These 
calculations  enhance  physical  understanding  and  help  to  dispel  con- 
fusion about  this  reflection  process.  In  particular,  it  is  shown  that 
the  mechanism  which  ionizes  the  gas  plays  no  role  in  the  wave  energetics, 
and  that  (unlike  the  case  of  scattering  off  a moving  plasma)  no  energy 
is  exchanged  between  the  wave  and  the  plasma  recoil.  Wave  energy  is, 
however,  deposited  in  the  kinetic  and  magnetic  energy  associated  with 
stationary  transverse  electron  currents  that  are  set  up  in  the  plasma 
as  a result  of  the  reflection  process. 

In  Section  IIC,  the  case  of  an  oncoming  front  moving  faster  than 
the  speed  of  light  is  treated.  Such  a front  is  physically  realizable, 
but  paradoxes  are  avoided,  as  it  is  shown  that  no  wave  reflection  occurs. 

In  Section  III  the  reflection  problem  is  formulated  and  solved  with 
no  restriction  on  the  shape  or  width  of  the  electron  density  profile  in 
the  ionization  front  (but  still  for  the  case  of  normal  incidence  on  an 
infinite  planar  front).  It  is  shown  that,  in  the  overdense  case, 

4 


W 
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Eqs.  (l)-(4)  for  U < 0,  and  (l),  (2),  (5),  (6)  for  U > 0,  hold  quite 
generally.  In  the  underdense  case,  explicit  solutions  for  the 
reflection  coefficient  are  written  in  terras  of  the  reflection  co- 
efficient for  a stationary  electron  density  profile  of  the  same  shape, 
in  a stationary  plasma.  In  Section  IIIB,  a complete  and  explicit 
solution  is  found  for  a particular  electron  density  profile  in  the 
front  (hyperbolic  tangent),  whose  width  can  be  varied  from  zero  to 
infinity.  In  Section  IIIC,  energy  conservation  is  demonstrated  and 
discussed  for  the  limiting  case  of  a broad  front,  based  on  a WKB 
calculation.7 

In  Section  IV,  the  theory  is  extended  to  electromagnetic  waves 
propagating  in  a waveguide,  and  to  waves  obliquely  incident  on  the 
front  in  free  space. 

A considerable  theoretical  literature  exists  on  reflection  at 
moving  discontinuities  in  material  properties, a~15  most  of  which  is 
incorrect  when  applied  to  plasma  ionization  fronts.  In  Section  V, 
we  briefly  comment  on  this  literature,  and  discuss  in  some  detail 
the  macroscopic  model  of  electromagnetic  scattering  at  a moving  dis- 
continuity in  the  dielectric  coefficient  of  a stationary  dielectric.12 
It  is  shown  that  (contrary  to  some  claims  in  the  literature13^14) 
this  model  is  inapplicable  to  moving  ionization  fronts,  and  is,  in 
fact,  only  one  of  several  possible  models  for  a discontinuity  in  a 
dielectric.  These  models,  which  lead  to  different  reflection  co- 
efficients, must  be  sorted  out  on  the  basis  of  the  microscopic  physics. 

In  Section  VI,  we  discuss  the  technological  applications  of  the 
ionization  front  scattering  process. 
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Section  VII  summarizes  our  results. 


II.  Sharp  Ionization  or  Recombination  Front 

A.  Solution  of  the  Scattering  Problem 

We  shall  perform  our  calculation  in  the  frame  of  reference  in 
which  the  ionization  or  recombination  front  is  stationary  (the  "front 
frame").2  In  this  frame,  in  the  ionization  front  case,  neutral  gas 
flows  into  the  discontinuity  from  the  left  at  velocity  U > 0,  and 
plasma  flows  away  to  the  right  at  the  same-6  velocity  U,  as  illustrated 
in  Fig.  1.  In  the  case  of  a recombination  front,  the  only  difference 
is  that  U <0,  i.e.  plasma  flows  into  the  front  and  gas  flows  out. 

In  the  front  frame,  the  scattering  problem  is  stationary  and  thus  an 
incident  wave  of  frequency  u>  generates  reflected  and  transmitted  waves 
at  the  same  frequency.  We  shall  assume  throughout  this  paper,  for 
convenience,  that  the  neutral  gas  has  the  dielectric  properties  of 
the  vacuum.  In  the  present  section,  we  assume  that  the  front  is  much 
narrower  than  the  wavelength  of  any  wave,  so  that  the  electron  density 
n(z)  in  Fig.  1 can  be  represented  as  a step  function,  discontinuous 
at  the  front.  In  Section  III  we  shall  consider  the  general  case  where 
n(z)  is  arbitrary. 

We  first  solve  for  the  transverse  normal  modes  in  the  plasma, 
where  the  electron  density  n is  constant.  The  transverse  electron 
velocity  v in  the  presence  of  the  transverse  wave  fields  E and  B is 
determined  by  the  linearized  momentum  conservation  equation  i.e. 
generalized  Ohm's  law), 


b 


+ u 


(7) 


5v 

TT 


fi v 
Jz 


— + — U x b\ 

V V C ~ 7 


In  Eq.  ( 7 ),  Y = (1  - @2)-1^2.  Also,  using  Maxwell's  equations  and 
neglecting  the  perturbed  ion  velocity. 


fiB 


V x E = - - 


V x B = - 


c fit  ’ 
SE 


(8) 


l4.tr 


nev 


c fit  c 

and  assuming  spatial  and  temporal  dependence  exp(ikz  - itrt),  Eqs. 
(Y)_(9)  lead  to  the  dispersion  relation 


(9) 


(10) 


where  uu p2  - It-  nne2/mey.  Thus  there  are  three  normal  modes  at  any  given 


frequency  m. 


ki(m)  = 


CJL) 

C 


/ 0)  ; 

1 - -4 

\ U ) 


2 \l/2 

H > 

if  m2  > cu  2 

1/2 

P 

1 , 

if  a)2  < a)  2 , 

(Ha) 


k2  (oj)  = - kj.(<n). 


k (<u)  = (u/U 
m 


(lib) 

(11c) 


The  roots  kj.  and  k2  are  just  the  usual  rightward-  and  leftward- 
propagating  electromagnetic  waves  if  the  plasma  is  underdense,  i.e. 


(jUp  < m,  or  the  equally  well-known  evanescent  and  growing  waves  if  the 
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plasma  is  overdense,  i.e.  > u).  We  note  that  the  dispersion 
relation  for  electromagnetic  waves  is  covariant  in  the  form 


0)2  = (JD2  - k2c2. 


It  follows  that  cu  is  a scalar  under  Lorentz  transformation,  and 
P 

thus  is  numerically  equal  to  its  value  in  the  laboratory  frame, 

'ju*  = (krrn*e2/m  z . Hence,  in  the  laboratory  frame,  the  condition 
for  the  plasma  to  be  overdense  is 


•s  > •?  ten 1/2 


The  other  normal  mode,  k , is  often  overlooked  in  discussions 

of  transverse  waves  in  plasmas.  The  nature  of  this  wave  is  elucidated 

by  noting  that  E = (cu/kc)B  = 5B  , so  that  if  we  Lorentz  transform 

to  the  laboratory  frame,  where  the  plasma  is  stationary,  we  find 

E*  * a)*  =0.  Thus  in  the  laboratory  frame  this  "wave"  is  just  a cold 
m m 

z » A z 

stationary  magnetic  disturbance  with  B*  = B e and  J*  = ik*zxB  e , 
which  is  indeed  a solution  of  Maxwell's  equations  and  the  linearized 
fluid  equation  of  motion,  and  is  therefore  a permissible  mode  in  a cold 
collisionless  plasma.  In  a frame  where  the  plasma  is  streaming,  this 
mode  becomes  a genuine  undamped,  dispersionless  electromagnetic  wave, 
carrying  energy  and  momentum  and  propagating  at  the  plasma  streaming 
velocity. 

The  following  argument  illustrates  why  the  magnetic  wave  is 
excited  at  the  sharp  ionization  front.  An  electron  moving  in  the 
plasma  region,  in  the  presence  of  the  evanescent  wave  (we  assume, 


for  the  moment,  that  the  plasma  is  overdense),  obeys  the  equation  of 


motion 


QV  6 

_ = ^(£  + |lIx£)  = = °>  c 0)(1  - |f)ei(kz-(Ut).  (14) 

Integrating  back  over  an  electron's  orbit,  and  imposing  the  initial 
condition  that  the  electron  is  born  (ionized),  at  z = z'and  time 
t-(z-z*)/U,  with  zero  velocity  in  the  laboratory  frame, 

v[z',  t-(z-z')/u]  = 0,  (15) 


we  find 


ieE (0,0 ) i ' z z"  \ 

v(z,t)  = ■ - jexpi(kz-u)t)  - exp[-iu>(t-  •S^“)]j  • 


Since,  in  the  case  of  a sharp  ionization  front  at  z =0,  all  electrons 
are  bora  at  z'  = 0,  the  current  J is  given  by  J = -nev,  with  v given 
in  Eq . (1 6).  The  first  term  in  Eq.  (1 6)  leads  to  the  current  associated 
with  the  evanescent  wave,  while  the  second  term  is  the  source  of  the 
magnetic  wave,  and  is  associated  with  the  boundary  condition  of  a sharp 
ionization  front. 

We  now  set  about  matching  waves  at  the  vacuum-plasma  interface, 
to  determine  the  reflection  and  transmission  coefficients.  We  shall 
see  that  the  physics  is  different  for  an  advancing  ionization  front 
(U  > 0)  and  for  a retreating  recombination  front.  We  treat  the  case 
U > 0 first. 


Equations  (11)  indicate  that  there  are  three  possible  waves  in 
the  plasma  at  frequency  tu,  but  the  wave  k2  is  ruled  out  by  its  unphysical 
behavior  at  x = + »,  since  it  is  growing  if  the  plasma  is  overdense,  or 
leftward-propagating  if  the  plasma  is  underdense.  Thus  at  the  interface, 
we  must  match  four  waves,  the  incoming  wave  (of  known  amplitude),  the 
reflected  wave  in  the  gas,  and  two  plasma  waves,  ki  and  k^.  Thus  there 
are  three  unknown  amplitudes  (one  more  than  in  the  usual  dielectric 
interface  problem),  and  we  need  three  matching  conditions.  It  is 
immediately  clear  from  integrating  Maxwell's  equations  (3,9)  across 
the  interface,  that  continuity  of  E and  B are  required,  as  usual.  A 
third  condition  is  supplied  by  noting  that  electrons  are  born  (by 
ionization)  at  the  interface  with  transverse  velocity  v = 0,  and  then 
stream  away  into  the  plasma  at  axial  velocity  U,  while  they  begin  to 
oscillate  in  the  perpendicular  plane.  Thus  the  current  at  the  inter- 
face is 


jJ(z  ■ 0)  =*  -nev(z  = 0)  =*  0, 

so  that  is  continuous  across  the  interface.  Equation  (9)  thus 
5B 

indicates  that  — r — is  also  continuous  across  the  interface. 

Oz 

By  noting  that 


(17) 


and 


B = (kc/<u)E, 


($B/5z)  = (ik2c/u;)E, 


( 15a ) 

(18b) 


using  Eqs . (11)  for  the  plasma  waves  and  = a)/c,  kr  = -au/c  for  the 
incident  and  reflected  vacuum  waves,  we  can  write  the  matching 
conditions, 


+ E + Ei, 

(19a) 

+ 6"1  Em  + WEi, 

(19b) 

+ 6"2  E + W2Ei, 
m 

(19c) 

(1  - a Jp2/^2)1/2, 

if  m2  > uo  2 

P 

(20a) 

i(u)p2/m2-l  )x/2, 

if  CD2  < 0)  2. 

P 

(20b) 

The  solution  of  Eqs.  (19)  for  the  amplitudes  of  the  reflected  and 
magnetic  waves  are 


Er  = 1-6  1-W 

E . 1+6  1-W 

1 


(21a) 


2P2(1-W) 

( 1+6  ) ( 1-6W ) 0 


(21b) 


The  magnetic  field  amplitudes  of  each  of  the  waves  are  determined  by 
(21)  and  (lSa)  and  we  can  then  write  down  the  electromagnetic  energy 
flux  (i.e.  Poynting  power)  P,  for  each  wave: 


46 4 1-W 

( 1+6  )“  1-6W  * 


'22b) 
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is  the  (amplitude)  reflection  coefficient  for  a stationary  sharp 


ionization  front.  For  the  overdense  case,|r|=  1.  We  also  note,  for 
use  in  the  next  section,  the  zz  components  of  the  momentum  flux, 

F = - < E2  + B^>/Stt,  for  each  wave  (where  the  brackets  indicate  a 
time  average ) : 

Fr  = Pr/c,  (24a) 


For  simplicity,  (24b)  is  given  only  for  the  overdense  case  (i.e.  W 
imaginary ) . 

We  now  transform  these  results  back  to  the  laboratory  frame. 

The  Doppler  shift  in  going  from  the  front  frame  to  the  laboratory 
frame  is 

tv*  = ouy  ( 1+3  ) , 

r (25) 

UJ*  « (JUV  ( 1-3  ) , 


so  that 


r 


Since  the  number  of  photons  in  the  incident  wave  packet,  or  in  the 
reflected  wave  packet,  is  frame-independent,  the  wave  packet  energy 
transforms  as 


®i  = eiY^1-9 


®r  = erY('1+e^ 


(27a) 


(27b) 


Combining  Eqs.  (22a)  and  (27)  gives 


e* 

r 


!U* 


1-6 

1+6 


(23) 


for  the  energy  reflection  coefficient  in  the  laboratory  frame.  Finally, 
the  pulse  duration  transforms  as  the  inverse  of  a frequency,17 


r 

T 


(29) 


whence  the  power  reflection  coefficient  in  the  laboratory  frame  is 
seen  to  be 


(30) 


For  the  overdense  case,  T = 1,  and 

We  now  turn  our  attention  to  the  case  U <0,  l.e.  a retreating 
recombination  front.  In  this  case  the  magnetic  wave  is  an  incoming 
traveling  wave  from  + ®,  and  thus  is  ruled  out  by  causality.  Thus 
the  wave  k (evanescent  if  the  plasma  is  overdense,  rightward-propagating 
if  overdense  1 is  the  only  permissible  wave  in  the  plasma,  and  the  usual 
matching  conditions  of  continuity  of  E and  B at  the  interface  still 
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hold,  but  the  third  condition,  continuity  of  , does  not  hold 

for  this  case.  It  derived  from  the  fact  that  J(z  =0)  * 0 for  an 
ionization  front  at  z = 0.  However  for  a recombination  front,  electrons 
dying  (by  recombination)  at  z « 0 carry  current  at  the  instant  of 
recombination,  and  thus  J(z  *0)^0. 

Thus  the  matching  conditions  become  the  same  as  in  the  elementary 
problem  of  a stationary  interface, 


Ei  ■ -Er  + Ei' 

(31a) 

E.  * +Er  + iWEi, 

(31b) 

whence  the  coefficients  of  reflection  and  transmission  are  found  to  be 

P 

-T--  |r|a  , 

ri 

(32a) 

- 1-! rj 2 . 

l 

(32b) 

Transforming  Eqs.  (32)  back  to  the  laboratory  frame. 

exactly  as  in 

the  derivation  of  Eqs.  (28)  and  (30),  we  find 
e*  x*  , la, 

£_  = I i ri 2 - | r! 2 

^ x*  i+!  e i 1 ! ' 

(33) 

P*  (D*  2 

“pf"  “ |r|2 

(3M 

for  the  coefficients  of  pulse  energy  and  power  reflection  in  the 
laboratory  frame. 
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B. 


Energy  Conservation 

We  shall  now  demonstrate  that  our  solution  conserves  energy. 

We  continue  to  work  in  the  front  frame  of  reference.  We  note  first 
that  we  have  assumed,  in  the  linear  electrodynamics  calculation  of 
Section  II  A,  that  the  streaming  velocity  U is  constant,  and  equal 
in  the  plasma  and  gas.  In  fact  U does  change  across  the  front, 
since  the  plasma  must  recoil  to  conserve  momentum,  but  it  is  consistent 
to  neglect  the  recoil  velocity  in  a linear  calculation,  since  it  is 
second  order  in  the  wave  amplitude  (and  also  is  inversely  proportional 
to  the  ion  mass).  However,  the  recoil  must  be  included  to  satisfy  the 
conservation  laws  (in  the  front  frame).  The  situation  is  analogous  to 
bouncing  a light  ball  off  a very  heavy  moving  object  (e.g.  the  earth): 
the  recoil  velocity  of  the  earth  is  negligible,  yet  it  absorbs  the 
momentum  change  of  the  ball. 

The  recoil  velocity  AU  could  be  calculated  from  a second  order 
electrodynamic  treatment  (including  ponderomotive  force).  We  take 
the  simpler  point  of  view  here  that  AU  is  determined  by  momentum 
conservation;  we  can  then  go  on  to  verify  that  energy  is  also  con- 
served. We  also  note  that  a small  density  change  at  the  front  is 
associated  with  the  recoil,  through  the  continuity  equation.  However, 
quasineutrality  requires  that  the  electron  and  ion  densities  be  locally 
equal  (assuming  singly  charged  ions,  for  simplicity). 

With  this  preamble,  we  proceed  to  write  down  the  conservation 
equations  for  particles,  the  z-component  of  momentum,  and  energy, 
considering  first  the  case  U > 0 (ionization  front): 
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0 = A ( nU ) , 

(35a) 

Pi+Pr-Pm  = 

(35b) 

cP.-cPr-Fm  = A[nUc2(meYe+miYi)] • 

(35c) 

The  terms  on  the  left  hand  side  of  Eqs.  (35b, c)  are  the  electromagnetic 
momentum  and  energy  flux,  P and  F respectively,  for  the  incident, 
reflected,  and  magnetic  waves,  given  by  Eqs.  (22,24).  We  consider 
here  only  the  overdense  case,  in  which  no  energy  or  momentum  flux  is 
associated  with  the  evanescent  wave  in  the  plasma.  On  the  right 
hand  side  are  the  changes  in  particle  momentum  and  energy  flux  due 
to  the  presence  of  the  waves;  A can  be  interpreted  as  the  difference 
in  a quantity  in  the  plasma  from  what  it  would  be  if  there  were  no 
incident  wave.  The  right  hand  side  includes  the  changes  due  to  both 
the  recoil  velocity  AU,  and  the  transverse  electron  oscillations 
associated  with  the  magnetic  wave.  Thus 


= Y3UAU/c2,  (36a) 

Aye  = Y3(UAU  + i < v2  >)/c2,  (36b) 

where  y - (l-U2/c2 )-1/2  is  the  unperturbed  value  (same  for  electrons 
and  ions)  and  < v2  > is  the  mean  square  oscillating  velocity. 
According  to  Eqs.  (9)  and  (11a),  v is  given  by 


v 


(-32+l)E. 


(37) 
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Using  the  magnetic  wave  amplitude  (21b)  in  (37)  then  gives 


nm  v 
e 


< v2  > 


4(1-6  )2 

1-6  ^+8  “to  ^/u)^ 

P 


(38) 


Solving  Eqs.  (35a, b)  for  AU,  and  substituting  in  Eq.  (35c)>  verifies 
(after  some  tedious  algebra)  that  Eq.  (35c)  is  satisfied  identically, 
and  thus  energy  is  conserved. 

Three  points  concerning  this  energy  conservation  demonstration 
deserve  comment.  First,  the  external  mechanism  which  causes  the 
plasma  ionization  plays  no  role,  neither  contributing  energy  to  the 
waves  nor  absorbing  energy.  Consistent  models  could  in  fact  be  con- 
structed in  which  the  energy  associated  with  ionization  becomes 
negligibly  small.  Secondly,  part  of  the  incident  wave  energy  penetrates 
the  plasma,  in  the  form  of  electromagnetic  and  transverse  electron 
kinetic  energy  associated  with  the  magnetic  wave.  This  occurs  even 
though  the  plasma  is  overdense.  However,  the  situation  would  look 
rather  different  in  the  laboratory  frame,  where  the  plasma  is 
stationary  but  the  ionization  front  is  moving.  Since  the  magnetic 
"wave"  is  stationary  in  this  frame,  one  could  more  accurately  say 
that  energy  is  deposited  in  the  plasma  in  the  vicinity  of  the  ionization 
front,  and  remains  in  the  same  place  in  the  plasma  as  the  front  moves 
on.  Thirdly,  we  note  that  energy  transfer  to  the  plasma,  through  the 
recoil,  occurs  in  the  front  frame.  As  seen  in  the  laboratory  frame 
where  the  plasma  is  stationary,  the  plasma  kinetic  energy  of  recoil 
would  be  vanishingly  small  [ order (AU2)] . The  comparison  with  the  case 
of  a "moving  mirror",  i.e.  a physically  moving  plasma,  is  instructive: 
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here  the  recoiling  plasma,  which  is  moving  in  the  laboratory  frame, 
does  exchange  energy  with  the  wave  in  the  laboratory  frame;  hence 
the  possibility  that  reflected  wave  energy  exceeds  incident  wave  energy 
in  this  case. 

We  shall  now  briefly  discuss  energy  conservation  for  the  recom- 
bination front  case,  U <0.  The  formal  framework  of  the  calculation 
is  similar  to  Eqs.  (35)*  but  it  is  now  appropriate  to  consider  A to 
represent  the  changes  in  quantities  at  the  recombination  front,  due 
to  the  presence  of  the  waves.  In  this  case,  the  magnetic  wave  is 
absent.  However,  the  electrons  acquire  a transverse  velocity  v,  by 
passing  through  the  evanescent  wave  en  route  to  the  recombination 
front.  A calculation  similar  to  that  of  Eq.  (33)  shows  that  just  to 
the  right  of  the  ionization  front 

i nmgY  < v2  > = E^2/4tt  . (39) 

Completing  the  calculation, just  as  in  the  case  U > 0,  then  verifies 
that  energy  is  conserved.  The  reflected  energy  in  the  laboratory 
frame  is  less  than  the  incident  energy;  the  difference  goes  into  the 
electron  transverse  kinetic  energy* Eq.  (39);  upon  recombination,  this 
energy  either  dissipates  by  heating  the  gas  or  goes  off  as  recom- 
bination radiation,  depending  on  whether  recombination  is  through 
three-body  or  radiactive  processes . 

C.  Super luminous  Velocity  of  the  Front 

It  is  quite  possible  for  an  ionization  or  recombination  front 
to  move  across  a gas  at  or  above  the  speed  of  light.  For  example, 


f 


IS 


the  ionization  front  could  be  produced  by  a laser  or  electron  beam 
sweeping  across  the  gas  (much  as  the  trace  in  a cathode  ray  tube  can 
move  across  the  screen  faster  than  c).  It  seems  clear  physically 
that  reflection  cannot  occur  at  such  a front,  since,  in  the  case  of 
a receding  recombination  front,  the  incident  electromagnetic  wave 
pulse  cannot  catch  up  with  the  front,  and  in  the  case  of  an  oncoming 
ionization  front,  a reflected  pulse  would  be  immediately  caught  by 
the  front.  What  actually  happens,  in  the  latter  case,  is  that  there 
is  no  reflected  wave,  but  two  propagating  waves  are  generated  in  the 
plasma,  as  will  now  be  shown. 

Since  U > c in  this  case,  we  shall  work  in  the  laboratory  frame. 
The  three  normal  modes  in  the  stationary  plasma  are  then  the  magnetic 


wave  (cjj#  * 0 ),  the  rightward-propagating  or  evanescent  (depending 
m 


on  whether  wf  > w^)  wave. 


k*  (ui£)  * (uj*/c  ) ( l-‘jjp2/x*2  )x/2 

and  the  leftward-propagating  or  growing  wave. 


(i*Oa) 


fcS(iD2)  = -(ui^/c)(l-(vp2/a^2)1^2.  (tob) 

Since  the  reflection  problem  is  non- stationary  in  the  laboratory 
frame,  both  k*  and  ji*  for  the  transmitted  waves  can  differ  from  the 
incident  wave  values.  However,  the  continuity  conditions  at  the  front 
can  be  met  only  if  all  waves  are  in  phase  at  the  front  (z  = -Ut),  i.e. 
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-k*Ut  - m*t  = -k*Ut  - w*t 
i i 


for  each  wave.  [Applying  this  condition  to  the  reflected  wave  is 
an  alternate  method  of  deriving  the  double  Doppler  shift  for  the 
sub luminous  case,  Eq.  (1).]  Thus 


-«l*(3+1)  = -k*U 
i m 


/ U)2  .1/2 

= -«£  ^1  - +1 

r / m2  x1/2  I 
= <4  Sy1  ' - ij  > 


where  0 = U/c  > 1.  Equations  (42b, c)  can  be  solved  for  uo*  and  o)^, 


with  the  result 


“?(«) 


_L_  + _L. 

3-1  + 3-1 


2 -.1/2 


Thus  auf  and  m*  are  both  real,  which  implies  that  kf  and  k^  are  also 

real  [from  the  phase-matching  conditions  (42b, c)];  it  follows  that 

neither  wave  is  cut  off,  l.e.  id  2 < m*2  and  u ? < uo*2.  Moreover  one 

pi  p 2 

of  the  waves  is  rightward-propagating  and  the  other  is  leftward- 
propagating  (i.e.  runs  after  the  front,  but  never  catches  up  with 
it).  Thus  the  three  amplitude-matching  conditions  discussed  in 
Sec.  IIA  (continuity  of  E,  B,  and  &B/6z)  can  be  met  by  choosing 
the  amplitudes  of  the  three  permissable  plasma  waves,  with  no 
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reflected  wave  (we  omit  the  statement  of  the  results  here).  The 
interpretation  of  the  leftward  propagating  wave  is  that  it  is  excited 
at  the  boundary  and  is  then  left  behind  by  the  front.  It  does  not 
represent  energy  launched  by  a source  at  z = + <*>  and  so  the  usual 
Sommerfield  radiation  condition  does  not  apply. 

Ill . Continuous  Transition 

We  now  consider  an  ionization  or  recombination  front  of  arbitrary 
shape  and  width,  but  with  U < c and  normal  plane  wave  incidence  from 
the  gas.  We  work  in  the  front  frame,  wherein  the  electron  density 
n(z)  is  stationary,  and  the  plasma/gas  streaming  velocity  U is  constant 
in  space  and  time.  For  the  ionization  front  case  (U  >0),  the  mean 
transverse  velocity18  of  the  electron  fluid,  v(z,t),  is  determined  by 
the  momentum  conservation  equation. 


nmeY  +mevu  _ 


= nmeY»—  + 


&v  \ 

uinr) 


Ull  1 

+meYUv  = - ne  (E  + - UxB).  (44) 


Ionization  does  not  constitute  a source  of  electron  transverse  momentum, 
since  electrons  are  born  with  v = 0.  Thus,  as  new  electrons  are  born, 
they  share  in  the  local  transverse  electron  momentum,  thereby  tending 
to  reduce  y.  This  is  the  interpretation  of  the  term m^Uvdn/dz  on  the 
left-hand  side  of  Eq.  (44).  Using  (44)  and  Maxwell's  equations  (8,9), 
and  requiring  that  all  quantities  have  the  time  dependence  exp(-iu;t) 
of  the  incident  wave  (since  the  medium  is  stationary  in  the  front 


frame),  we  arrive  at  the  following  wave  equation, 
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[c2d2/dz2  + u)2  - o)2(z)](-ioj  + Ud/dz  )E  = 0, 


(45) 


The  case  of  a recombination  front  is  somewhat  different.  Since 
the  transverse  momentum  of  an'  electron  is  lost  from  the  electron  fluid 
when  it  recombines,  a momentum  sink  mvvp(z)  must  be  added  to  the  right- 
hand  side  of  the  momentum  equation,  where  P — Udn/dz  is  the  local  recom- 
bination rate.  Thus  the  equation  for  v becomes 


nm  v 

e 


et 


+ u 


Qv 

9z 


= -ne (E  + - U x B ) , 
~ c ~ ~ 


(46) 


Maxwell's  equations  and  (46)  lead  to  a different  wave  equation. 


-ix  + Ud/dz  )uT2  (c‘=d2/dz2  + cu2  - au  2)E  = 0, 
P P 


(47) 


for  the  recombination  front  case. 

The  recombination  front  case  is  most  easily  treated.  Equation 
(4-7)  can  be  immediately  integrated  once  to  yield 


(c 


d2 

2 2 + u)2  - u)  2 )E  = E (u  2 exp(ioiz/u). 


dzc 


o p 


(48) 


For  large  positive  z,  au^  becomes  constant  and  Eq.  (48)  possesses  the 
following  particular  solution  : 


E =*  oj  2 E [m2(l-8-2)-'jD  2]-1  exp  i(  ^ z),  z -»  + °°. 
p p o p u 


(49) 


However,  for  the  recombination  front  case,  U < 0 and  (49)  represents  a 
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magnetic  wave  transporting  energy  in  from  z = +°°.  Since  this  is 
physically  unacceptable,  the  constant  Eq  must  be  zero,  and  (48)  reduces 
to  the  usual  equation  for  wave  propagation  in  a stationary  plasma. 


Thus  the  problem  of  reflection  from  a receding  recombination  front  is 
equivalent,  in  all  respects,  to  reflection  from  a receding  physically 
moving  plasma  with  the  same  electron  density  profile  (in  the  frame  in 
which  the  profile  is  stationary).  In  the  recombination  front  frame, 
the  solution  to  (50),  in  the  gas  (z  ■*  -®,  oj^  -»  0),  is  of  the  form 

E(z)  = E(eiajz/c  + T e”ituz/c).  ' (51) 

F is  the  (amplitude)  reflection  coefficient  for  the  stationary  front 
in  stationary  plasma.  It  is  given  exactly  by  Eqs. (20 ) and  (25)  if  the 
front  happens  to  be  abrupt  (which  means  that  its  width  is  much  less 
than  a wavelength).  It  is  also  well  known  that  1 1”  j = 1 if  the  plasma 
is  overdense.  In  other  cases,  the  calculation  of  T is  a well-known, 
but  usually  complicated,  mathematical  problem.  In  general,  the 
reflected  power  in  the  front  frame  is  given  in  terms  of  T by  Eq.  (52a), 
and  the  reflected  power  and  pulse  energy  in  the  laboratory  frame  by 
Eqs.  (53)  and  (54-).  The  results  of  the  previous  section  now  become  a 
special  case,  for  which  T is  known  exactly,  of  this  general  result. 

We  can  also  write  down  the  general  solution  for  the  ionization 
front  case,  Eq . (45): 


— - . 
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E 


z 


(52) 


= j dz'  g(z')exp  ji(  ^ 


J 


where  g(z)  is  a solution  to  the  stationary  medium  reflection  problem 
(50)  with  appropriate  outgoing  wave  boundary  conditions,  and  e is  a 
positive  infinitesmal  introduced  for  convergence  (e  "*  0+).  In  writing 
(52),  a particular  solution  Ae*  z which  would  appear  on  the  right 
hand  side  has  been  neglected  (A  * 0),  since  A f 0 would  imply  the 
existence  of  a rightward-propagating  magnetic  wave  component, 

E ~ exp ( ituz/U ),  in  the  plasma- free  region  z -*  where  the  incoming 
wave  is  specified  to  be  an  electromagnetic  wave  of  given  amplitude, 

E = Ei  exp(iujz/c  ) . 

We  notice  that  to  find  the  reflection  coefficient,  we  need  to 
evaluate  Eq.  (52)  only  asymptotically  at  z •*  and  that  to  do  this 
we  only  need  the  solution  g(z)  of  (50)  asymptotically  at  z ■*  - ®, 

But  the  latter  is  given  by  Eq.  (51),  in  terms  of  the  constant  T (the 
reflection  coefficient  for  the  stationary  problem).  Substituting  (51) 
in  (52),  we  can  perform  the  integration  exactly,  and  find 


The  reflected  power,  in  the  front  frame,  is  then 


I 


It  is  thus  generally  true  that  the  reflection  coefficient  (in  the 
frame  of  a moving  ionization  front)  is  ( 1-3  '2/ ( 1+3  )2  times  the 
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reflection  coefficient  for  a stationary  plasma  with  the  same  density 
profile  (and  thus  is  always  less  than  unity).  The  results  of  Sec.  IIA 
are  a special  case  of  this  conclusion.  In  particular,  if  u)2  < u^2, 
then  |T|2  = 1 and  the  reflection  coefficient  is  independent  of  the 
exact  profile  n(z)! 

Lorentz  transforming  to  the  laboratory  frame,  and  following  the 
same  reasoning  as  in  Sec.  IIA,  we  find  the  reflected  power  and  pulse 
energy  to  be 


(55) 


|r|2  < i. 


(56) 


B . Exact  Solution  for  a Particular  Profile 

Having  obtained  a rather  general  result  for  the  magnitude  of  the 
reflection  coefficient,  we  now  consider  the  excitation  of  the  magnetic 
wave,  discussed  in  Sec.  II  for  the  sharp  boundary  case.  We  shall  find 
that  the  situation  is  radically  different  for  the  case  of  a sharp 
boundary  and  for  the  case  of  a broad  transition  many  wavelengths  thick. 
In  order  to  illustrate  this  we  consider  a specific  profile 


CD2  (z)  = i x2o  [l  + tanh  (z/2a)]. 


(57) 


and  assume  that  u>  2{<°)  = a^2  > w2,  i.e.  the  plasma  is  overdense.  The 
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exact  solution  g(z)  to  Eq.  (50);  with  this  functional  form  for 
is  known  in  analytic  form,*9  in  terms  of  the  hypergeometric  function. 

By  using  the  Barnes  integral  representation  of  the  hypergeometric  function, 
the  expression  for  g(z)  in  Ref.  19  may  be  written  in  the  following  form1": 


g(z) 


-akz 

e 

2TTi 


I 


ds 


r (Ai  )r  (Ag  )r  ( 1-3  )exp[  z(  l-s  )/a] 


F(s  + 2kra  ) 


(53) 


where  Ai ; a = (s-1  + ikr  + kcar/,3- 


(a,2 

po 


- u)2)1/2/^,  r(w) 


are  gamma 


functions,  g (z)  ~ exp. / -akz)  for  large  z/r  as  required,  and  the  contour  C 
is  illustrated  in  Fig.  2.  Since  F(w)  has  no  zeros  and  an  infinite  number 
of  poles  at  w * 0,  -1,  -2,...,  with  residues  (-l)w/(-w) I,  the  integrand 
of  Eq.  (50)  has  the  poles  labeled  with  x's  in  the  figure.  The  series 
of  poles  on  Im  s)  * kc  and  on  Im( 3 ) = - kc  result  from  I*(Ai)  and  F (A2 ) 
respectively;  while  the  series  of  poles  on  the  real  axis  result  from 
F l-s).  Putting  (5®)  in  (52)  and  interchanging  the  orders  of  the  z' 
and  s integrations, we  obtain 


ce 


-ksrz 


2*ri 


/ 


ds 


"(Ai  )F (As )F ( l-s )exp[ - ( s-1 )z/ct] 
F ' s+2k o or  " s - 1+k  ta  + ixc/U  * 
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Thus  another  pole  is  introduced  in  the  integrand  at  s 3 l-kffcr  - ixc/u 
(labeled  by  a dot  in  Fig.  2).  Note  that  the  contour  C lies  to  the 
left  of  this  pole  so  that  (as  we  shall  see)  E satisfies  the  condition 
of  having  no  component  proportional  to  exp ( ixz/U ) in  the  region 
z « - c.  For  z > 0,  the  contour  can  be  closed  with  a large  semicircle 
at  infinity  in  Re(s)  >0.  For  z » a the  dominant  contribution  is  from 
the  pole  at  s * 1 -kc  y + ixc/U,  with  all  other  pole  contributions  going 
exponentially  to  zero.  Thus 
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E — X3  exp(iu«/U),  for  z » c, 


60a ) 


where  X3  is  obtained  from  the  residues  at  the  pole.  Similarly  we 
obtain 

E 2=  X!elkz  + X2e‘ikz,  for  z « -a  , (60b) 

which  results  from  the  dominant  pole  contributions  at  s = 1-ko  a + ik a . 
Evaluating  the  residues,  Xi,  X2  and  X3,  we  obtain  the  reflection  co- 
efficient | X2/Xx l 2 = ( 1-3  )2/( 1+6  )2,  in  accordance  with  the  general  result 
(54).  We  also  find  the  excitation  coefficient  for  the  magnetic  wave 

2tt32  1-3  ( Tpo\2  2k o sinh(2TTk o ) 

S^a^+l  1+6  \ u)  / sinhlfTka  ( 1-3  )/S]  sinhlTTkc ; 1+3  )/s!]  * 

(61) 

In  obtaining  (6l)  we  have  made  use  of  the  following  gamma  function 
identities:  |F(iy)!2  = n[y  sinh(TTy)]-1  ;where  y is  real),  T (z+l)  = 
zf(z),  and  j T ( 2 ) | 2 = |rfz*)|2  (where  z * is  the  conjugate  of  z).  For 
a -*  0,  Eq.  (6l)  yields  the  result  for  the  sharp  boundary  case  found 
in  Sec.  II.  For  a -*  «,  the  excitation  coefficient  for  the  magnetic 
wave  becomes  exponentially  small.  While  we  have  demonstrated  this 
only  for  a particular  profile,  this  same  result  can  be  shown  in  a 
general  way  for  a large  class  of  profiles.  This  may  be  interpreted 
as  follows:  When  a particle  is  suddenly  ionized,  its  transverse  velocity 

acquires  an  oscillatory  component  at  the  local  frequency  of  the  usual 
electromagnetic  wave,  and  also  a component  which  has  zero  frequency  in 


X3 

xT 
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the  laboratory  frame  (as  discussed  in  Sec.  II ).  For  the  case  of  a 
sharp  front,  this  transverse  streaming  part  created  at  the  front  sets 
up  currents  which  support  the  magnetic  wave.  For  the  broad  front, 
however,  the  transverse  streaming  velocity  component  of  electrons  is 
created  continuously  throughout  the  broad  front.  At  any  point  behind 
the  front  there  will  be  electrons  bom  at  different  points  in  the  front, 
having  correspondingly  different  transverse  velocities.  The  velocities 
of  these  electrons  add  with  different  phases  to  produce  a reduced  current, 
and  when  the  front  is  very  broad,  this  phase  mixing  results  in  an 
exponentially  small  magnetic  wave. 

The  analyses  of  Sec.  II  and  of  the  present  section  present  us  with 
the  following  puzzle.  For  the  sharp  boundary  case  we  found  (Sec.  II) 
that  the  energy  lost  upon  reflection  was  used  to  excite  the  magnetic  wave. 
However,  for  the  broad  boundary  case  we’ have  now  shown  that  the  excitation 
of  the  magnetic  wave  becomes  exponentially  small.  Yet  we  have  shown  that 
the  reflection  coefficient  is  the  same  in  both  cases.  In  the  next  sub- 
section we  resolve  this  apparent  inconsistency. 

C.  Energy  Conservation  for  the  Case  of  a Broad  Ionization  Front 
(c  > U > 0) 

As  before  we  work  in  the  front  frame  and  attempt  to  demonstrate 
conservation  of  energy.  The  key  point  for  the  limit  of  a broad  front 
is  that  even  though  there  is  no  transmitted  magnetic  wave,  energy  is 
still  transferred  into  the  plasma  by  the  phase  mixed  transverse  stream- 
ing electron  velocities-  i.e.  a transverse  electron  temperature  is 
produced.  As  seen  in  the  front  frame,  energy  transfer  also  occurs 
through  the  plasma  axial  recoil  velocity  !U . 
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The  conservation  equations  for  axial  momentum  and  energy  are 
similar  to  those  discussed  in  the  sharp  front  case , Eqs.  (35 )>  except 
for  the  absence  of  the  electromagnetic  fluxes  associated  with  the 
magnetic  wave,  i.e. 


Pi+Pr  “ ^CnU2(meYe+m1Yi)] 

= i ra  nv3@2  < v2  > + (m  +m. )nY3UAU, 

2 e e 1 ’ 

cPt-cPr  * A[nUc2(meYe-HniVi )] 

= i n^nUY3  < V2  > + (me+m^ )nY3U2AU, 


(62a) 


'62b) 


where  A again  denotes  the  difference  between  the  value  of  a quantity 
in  the  region  z •*  + <*>  with  and  without  the  presence  of  the  incident 
electromagnetic  wave  [note  that  A(nU)  = 0] . Eliminating  AU  from  (62) 
and  using  (54),  we  obtain  the  condition  for  energy  conservation. 


cPi  = k "V1'*'3  ( 1+9  )2  c < v2  > . 


(63) 


It  thus  remains  to  verify  that  this  energy  conservation  condition  is 
indeed  satisfied  for  the  case  of  a broad  front.  To  calculate  g;'z)> 
the  solution  of  (50),  for  a very  gradual  transition,  it  is  adequate 
to  use  the  WKB  approximation7  for  g(z), 

g(z)  s=  Ad-m  2(z )/o)2 )”1/4[ei^  + e-1*  + ®] , if  m2  > 2(z),  (64) 


g(z)  - 0,  if  xs  < uj2(z  ), 


(65) 
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where  'Jr  = ik  J | 1_lu)  (z  * )/uj]2  | dz'  and  9 is  a constant  phase 
factor.  Using  [6k)  and  (65)  in  (52),  we  have  ['ju2  > ou2(z)] 


A 

/ V/4 

1-  J- 


(1-0 ) eiv 

, US 

/ V 

8(1-  7^ 

1+8  Jl — -s^—J 

\ 5 / 

V x / 

e-ime  (66' 


The  transverse  velocity  of  an  electron  created  at  time  t and  position  z 
can  be  regarded  as  the  sum  of  a streaming  component  (zero  frequency 
in  the  laboratory  frame)  and  an  oscillating  component;  the  total 
electron  transverse  velocity  at  the  instant  of  ionization  must  be  zero. 
By  integrating  the  transverse  acceleration  due  to  both  the  incident  and 
reflected  waves,  we  can  calculate  the  transverse  streaming  velocity  com- 
ponent; we  then  find 


m nv3  < v2  > * m v3  | — -j  2 ( 1-32 ) f ° l 1-u > 2(  z )/u>2]  _1/2 
e e 1 mur/  1 P ( o 

{[  1-3  (l-u>2(z  Vu>2  )1/2]-2  + [ 1+9  ( 1-u)2 ' z Vuj2 "I1  ^2] ”2 1 dn( z )/dz )dz, 

' D P 


where  z is  defined  by  ou  2'z  ) * x2.  Making  a change  of  integration 
o p o 

variable  from  z to  W = ( l-x2(  z )/x2  )1/,;2,  we  readily  obtain 


. 2 

T^mgny^l+6^  c <v^> 


c 

Ltt  ’ 


which  is  the  same  as  cP^,  completing  the  proof. 
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IV.  Cases  of  Oblique  Wave  Incidence  and  Incidence  in  a Waveguide 

If  the  incident  wave  propagates  obliquely  to  an  infinite  planar 
ionization  front,  or  if  the  waves  and  plasma  are  confined  within  a 
waveguide,  the  waves  have  a perpendicular  component  k^  of  the  propa- 
gation vector.  In  this  section,  we  extend  the  results  of  the  previous 
section  by  including  kx  in  the  calculation.  We  consider  only  the  case 
where  g is  transverse  to  the  electron  density  gradient  Vn^,  i.e.  E 
is  perpendicular  to  the  plane  of  incidence  for  the  oblique  incidence 
case,  or  the  incident  mode  is  a transverse  electric  (TE ) mode  in  the 
waveguide  case.  The  other  polarization  leads  to  the  well-known 
complications  of  resonant  absorption  (e.g.  excitation  of  electrostatic 
Langmuir  waves  ) even  for  U = 0,  and  is  presently  under  study. 

Again  working  in  the  front  frame,  from  Maxwell's  curl  E equation, 
U x B/c  = i(r_1UdE/dz.  Also  from  Maxwell's  equations  we  have 

(d2/dz2  + jj~/cz  - k^2)E  = ^fTnefcv/bt.  (69) 

Here  k^  = (uu^/c)  sin  9/*  for  the  obliquely  incident  plane  wave  case, 
with  6j*  the  angle  of  incidence  (in  the  laboratory  frame).  For  the 
waveguide  case,  k^  is  the  usual  waveguide  mode  transverse  component, 
e.g.  k2  = (pn/a)2  + (qrr/b)2  for  the  pq  mode  of  a rectangular  a by  b 
waveguide.  Using  these  equations  together  with  (UU)  yields 

(-iu)  + Ud/dz)o)”2  (c2d2/dz2  + uj2  - cu  2 - k2c2  'E  = 0 (70 ) 

P P 

for  the  recombination  front  case,  and 
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(c2d2/dz2  + x2  - x2  - k 2c2)(-ix  + Ud/dz)E  = 0 
P x 


(7D 


for  the  ionization  front  case.  These  equations  differ  from  Eqs. 
(45^47)  only  by  the  inclusion  of  the  k2c2  term.  Following  the  rea- 
soning used  in  Sec.  Ill,  we  again  find  that  the  recombination  front 
reflection  coefficient  (in  the  front  frame)  is  the  same  as  for  a 
stationary  front  (U  =0),  i.e. 

ik  z -ik  z 

E(z  -*  -”)  = E^e  2 + Fe  2 ) (72) 

where  T is  again  the  reflection  coefficient  for  a stationary  front 
(U  =0).  For  the  ionization  front,  we  find,  similarly  to  Eq.  (53 ), 


ik  z 1-k  U/x  -ik  z 

E(z  -*-«)■  E.  (e  2 H . . |2'  > — r 2 ), 

' iv  * l+k2U/x  e 

where  x and  k2  are  the  frequency  and  z-component  of 
in  the  front  frame. 

Lorentz  transformation  to  the  laboratory  frame 
of  reflected  to  incident  frequencies  can  be  written 
equivalent  forms 


(73) 

the  wave  number 

shows  that  the  ratio 
in  any  of  the  three 


to  * 1+k  U/cd 

r __  2 

x^  1-k  U/x 
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(71**) 

(74b) 

r 
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TTp > 

(74c) 
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where  6*  and  9*  are  the  angles  of  incidence  and  reflection  in  the 
laboratory  frame.  Equations  (7^a)  and  (7^b)  exhibit  admirable  symmetry, 
but  only  (7kc)  is  completely  defined  in  terms  of  incident  wave  para- 
meters. Thus  Eqs.  (74)  can  also  be  used  to  solve  explicitly  for  k^ 
and  9*. 

It  follows  from  Eqs.  (72)- (74a)  that  the  ratio  of  incident  to 
reflected  power,  in  the  front  frame,  is 

Pr/Pt  - |r!2  (75) 

for  a recombination  front,  and 


Pr/Pt  = (x*/w*)a|r|2 


(7 6) 


for  an  ionization  front.  As  in  Secs.  II  and  III,  we  find  the  ratio 

of  reflected  wave  packet  energy  e*  to  incident  packet  energy  e*,  in 

the  laboratory  frame,  by  using  Eqs.  (74-76)  and  arguing  that  wave 

action,  i.e.  number  of  photons,  is  a Lorentz  invariant,  so  that 

e , /uj - = e*/(u*  and  e /a)  = e*/x*.  We  then  find 

ii  ii  rr  rr 


e * 
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lrl2 


'77 ) 


for  a recombination  front,  and 


(78) 
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for  an  ionization  front.  Equations  (75-78)  are  identical  to  the 

results  for  normally  incident  plane  waves,  and  the  Doppler  shift  (74) 

differs  only  in  the  substitution  of  k for  k. 

z 

All  of  these  results  apply  to  both  the  obliquely  incident  plane 
wave  on  an  infinite  planar  front,  and  to  incidence  in  a waveguide. 
However,  the  power  reflection  coefficient  in  the  laboratory  frame  is 
different  in  the  two  cases,  since 

Pr*/Pi*  = (er*/Tr*)/(ei*/ri*),  (79) 

and  the  pulse  length  ^depends  on  the  group  velocity,  which  differs 
in  the  two  cases.  For  the  case  of  oblique  incidence  in  free  space, 
the  group  velocity  equals  the  ph,  ’e  velocity, 


T_.*/Tr*  = U*/v*, 


and 


P * 
r 


i 


(81) 


just  as  in  Secs.  II  and  III.  However  for  a waveguide, 


V/Tr*  = (3gs-e)/(e+eg)> 


where' 


and 


9g  = (1  - ki2caAiK*'2)1/2 


1+8 2 + 2SS  “1 

9-  - 9„  ~ 


gs  g 1+3^  + 2SS 


g 


(82) 


(83) 


(84) 
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V . Critique  of  Previous  Work:  Moving  Interface  Between 

Dielectrics 

A considerable  (but  rather  disjointed)  theoretical  literature 
exists  on  reflection  and  transmission  of  electromagnetic  waves  at 
moving  discontinuities  in  electrodynamic  properties  of  media.8-15 
Although  many  of  the  authors  have  argued  that  their  calculations  apply 
to  moving  ionization  fronts,  the  existence  of  the  magnetic  wave  is 
overlooked  in  almost  all  of  these  works 12-15  Consequently  the 
boundary  value  problem  which  is  set  up,  and  the  resulting  solutions, 
are  incorrect  when  applied  to  plasma  ionization  fronts,  as  will  be 
discussed  below  in  some  detail.  We  note  that  the  sharp ly  discontinuous 
moving  ionization  front  problem  has  been  addressed  correctly  by  Hughes 
and  coworkersf1  However,  these  authors  were  interested  in  non- 
relativistic  shocks,  and  therefore  neglected  to  calculate  the  reflection 
coefficient,  which  differs  from  unity  only  by  terms  of  order  3 (and  con- 
stitutes our  principal  result).  Hughes,  et  al.^1  also  considered  the 
different  situation  of  a collision-dominated  plasma  with  collision 
frequency  large  compared  to  the  wave  and  plasma  frequencies. 

Of  the  remaining  calculations  in  the  literature,  several  set  up  the 
problem  of  reflection/transmission  at  a moving  discontinuity  in  dielectric 
coefficient.  Since  a strictly  dielectric  formulation  is  used,  it  is 
implicitly  assumed  that  the  medium,  like  a vacuum,  sustains  only  two 
transverse  modes  at  any  given  frequency,  and  the  possibility  of  a third 
mode  (i.e.  the  magnetic  mode  in  a plasma)  is  excluded.  We  shall  consider 
in  some  detail  the  model  of  Tsai  and  Auld,^2  wherein  the  moving  inter- 
face is  treated  as  a sharp  moving  discontinuity  between  stationary 


dispersionless  materials  with  different  dielectric  constant,  i.e.  it 
is  assumed  that,  within  each  region,  the  polarization  P(x, t)  is 
proportional  to  the  local  electric  field  E(x, t), 

P(x,t)  = x£(5,t),  (85) 

at  all  points.  Under  these  assumptions,  and  for  normal  incidence,  it 
is  found12  that  the  Doppler  shift  of  the  reflected  wave  is 


1+U/V  . 

ei-  , 

l-U/V  . 


the  pulse  length  is  reduced  as 

T * U)  * 

r i 

T * ' = ’ 

i r 


(86) 


(87) 


the  power  reflection  coefficient  is 


Pr*  ( \*  1-(®2/«i  )l/a  Y* 

TT  = l+(«a/«x)*>*/  ’ 

and  the  reflection  coefficient  for  total  pulse  energy  is 


(88) 


er*  “r*  / l-(*a/«i)l/ay 

V V \ l+(.a/«x)l/a/ 


(89) 


where  V ^ is  the  phase  velocity  of  the  incident  wave,  and  Si  and  S2 
are  the  dielectric  coefficients  on  the  two  sides  of  the  interface. 
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We  see  immediately  that  if  V ^ -*  c,  i.e.  the  dielectric 
susceptibility  is  weak  on  the  incident  side,  and  e2  is  associated 
with  the  plasma  dielectric,  l-u u^/x2,  Eqs.  (S6-69)  become  identical 
to  Eqs.  (1-4)  for  the  reflection  off  the  front  of  a moving  plasma.  In 
particular,  Eqs.  (86)  and  (89)  indicate  that  the  energy  of  the  reflected 
wave  pulse  may  exceed  that  of  the  incident  wave  pulse,  even  though  the 
dielectric  material  is  not  moving.  This  reflection  process  is  thus 
unlike  the  process  of  reflection  off  a moving  ionization  front  in  a 
stationary  plasma  for  which  we  have  seen  that  e^*  ^ e^*  always  holds. 

One  may  reasonably  wonder  about  the  reasons  for  this  difference,  and 
the  source  of  the  extra  reflected  energy  in  the  case  of  the  dielectric 
interface.  We  shall  not  perform  a complete  calculation  of  the  latter, 
but  we  hope  that  the  following  comments  will  help  to  clarify  these 
points . 

There  are  in  reality  no  completely  dispersionless  dielectrics  that 
always  obey  Eq.  (85).  Any  dielectric  may  be  considered  as  a collection 
of  oscillators,  with  natural  resonances  at  one  or  more  frequencies. 

Let  tu  (a  spatially  dependent  quantity)  be  the  lowest  such  resonance. 

Then  a number  of  tacit  assumptions  are  built  into  the  dielectric  dis- 
continuity model  of  Ref.  12: 

(a)  If  a is  the  width  of  the  transition  region,  then  c is  much 
less  than  the  incident  and  reflected  wavelengths,  since  a sharp 
transition  is  assumed. 

(b)  The  wave  frequency  x * at  any  given  point  is  much  less  than 
the  local  value  of  xq,  since  it  is  assumed  that  the  dispersionless 
relation  (85)  holds,  i.e.  that  the  dielectric  responds  to  the  wave 
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field  as  if  it  were  a time- independent  field. 

(c)  It  is  also  assumed  that 

<7  » U<u  "S  (90) 

since  this  is  the  requirement  for  the  change  in  dielectric  coefficient 
to  be  adiabatic,  so  that  Eq.  (35)  is  obeyed  near  the  interface.  Thus 
the  interface  cannot  be  a strict  discontinuity.  For  example,  let  us 
consider  the  simplest  possible  model  of  a dielectric:  a collection  of 

identical  classical  oscillators,  with  electrons  attached  to  atoms  by 
"springs".  Each  oscillator  then  obeys  the  equation 

x = - a:  2 x - eE/m  . i'9l) 

o e 

If  we  let  E be  constant  [^in  accordance  with  assumption  (a)  discussed 
above],  but  allow  ojq  to  vary  in  time,  as  the  front  passes  over  a given 
atom,  then 

x 81  - eE/m  id  2 (92) 

e o 

if  Eq.  (90)  holds;  this  then  allows  (35)  to  be  satisfied.  In  the 
opposite  limit  to  Eq.  (90 ),  e.g.  in  the  case  of  an  abrupt  transition 
from  u)  = 00  (tightly  bound  electrons  - zero  susceptibility)  to  finite 


33 


in  the  immediate  vicinity  of  the  interface,  i.e.  the  electrons  act  as 
if  they  were  free.  Thus  this  case  corresponds  microscopically  (within 
the  region  UuT1  from  the  front)  to  an  abrupt  ionization  front,  not  to 
a discontinuity  in  dispersionless  dielectrics. 

(d)  It  is  clear  from  these  arguments  that  for  the  model  of 
Ref.  12  to  hold,  each  individual  oscillator  must  make  a smooth  transi- 
tion in  accordance  with  Eq.  (90).  It  is  not  sufficient  for  the 
macroscopic  dielectric  coefficient  to  change  smoothly,  as  a result  of 
a statistical  distribution  of  many  oscillators,  each  of  which  makes  a 
sharp  transition. 

The  fact  that  the  reflected  energy  can  exceed  the  incident  wave 
energy,  according  to  the  dispersionless  dielectric  discontinuity  model 
of  Ref.  12,  can  be  explained  as  follows.  When  the  oscillator  fre- 
quencies are  changed  adiabatically,  in  presence  of  an  electric  field 
'due  to  the  waves),  work  is  done  on  the  oscillator  by  the  mechanism 
that  changes  its  frequency.  Thus  this  energy  is  available  to  enhance 
the  reflected  wave.  On  the  other  hand,  no  work  is  done  in  suddenly 
changing  the  oscillator  frequency.  These  results  can  be  established 
by  appropriately  integrating  Eq.  (91),  or  by  the  following  simple 
argument.  Imagine  an  electron  attached  by  a spring  to  an  atom,  in 
the  presence  of  an  electric  field.  If  the  spring  is  slowly  released, 
work  is  done  against  the  electric  field.  But  if  the  electron  is 
initially  tied  down  by  a string,  no  work  is  done  in  cutting  the  string, 
i.e.  instantaneously  increasing  the  susceptibility. 

Thus,  if  one  expects  to  apply  the  dispersionless  dielectric 
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discontinuity  model  to  a given  experiment,  one  must  verify  that 
assumptions  (a-d)  are  valid.  In  particular,  one  must  pay  attention 
to  assumption  (d),  which  deals  with  the  microscopic  physics. 

VI.  Applications 

We  have  shown  that  the  power  reflection  coefficient  for  upshifted 
microwaves  reflected  off  an  oncoming  overdense  ionization  front  is 
ideally  unity.  However,  the  reflected  pulse  is  shortened  as  in 
Eq.  (2);  the  total  reflected  wave  packet  energy,  e^*  = (uj^*/u)r*  )e^*, 
is  less  than  the  incident  energy  and  unfortunately  decreases  as  the 
Doppler  upshift  increases. 

Nevertheless,  upshifting  high-power  microwave  sources  by  factors 
of  10  to  100,  with  efficiency  of  10$  to  1$,  could  produce  very 
interesting  fluxes  in  the  millimeter  and  particularly  the  submilli- 
meter wavelength  regimes,  where  intense  sources  have  been  hard  to  come 
by.  This  is  particularly  true  if  a clean,  precisely  controllable 
moving  ionization  front  can  be  produced  without  an  energy  investment 
that  far  exceeds  the  energy  of  the  microwave  sources.  This  may  well 
be  possible,  since  the  ionization  energy  which  must  be  supplied  to 
form  an  overdense  plasma  is  quite  modest,  typically  of  the  order  of 
0.1  joule  per  meter  for  a 10cm  diameter  plasma  at  ng  - 1014cm-2. 

Many  techniques  for  producing  moving  ionization  fronts  may  be 
envisioned.  A high  current,  relativistic  electron  beam  propagating 
through  a gas  typically  creates  a moving  ionized  region  in  the 
stationary  gas,  with  electron  density  much  higher  than  that  of  the 
beam  itself.  It  is  also  possible  to  drive  a strong  electromagnetic 
pulse  down  a gas-filled  guide,  producing  an  ionization  front 
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accompanying  the  pulse.  Perhaps  the  most  promising  approach  for  our 
purposes  is  to  sweep  an  ionizing  laser  beam  across  the  gas  transversely. 
In  fact,  a variation  of  this  approach  has  recently  been  demonstrated 
in  connection  with  work  on  the  ionization  front  ion  accelerator.  In 
this  scheme,  a precisely  programmed  sequence  of  laser  pulses  is  fired 
(from  a series  of  many  light  pipes)  across  a cylinder  of  gas.  The 
laser  frequency  is  tuned  to  excite  the  gas,  which  is  then  ionized  by 
another  laser  propagating  axially  in  the  continuous  wave  mode. 

One  possible  limitation  on  the  microwave  power  that  can  be 
handled,  in  an  ionization  front  reflector,  would  arise  from  the 
requirement  that  the  incident  microwaves  themselves  do  not  pre-ionize 
the  gas.  Thus  a careful  selection  of  the  working  gas  and  ionization 
mechanism  may  be  necessary  to  optimize  such  a system. 

VII.  Summary  and  Conclusion 

In  this  paper  we  have  considered  the  problem  of  reflection  of  an 
incident  electromagnetic  wave  from  a relativistic  oncoming  ionization 
front  or  receding  recombination  front,  where  the  front  region  is 
assumed  to  separate  regions  of  stationary  neutral  gas  and  ionized 
plasma.  We  first  considered  the  case  of  a normally  incident  plane 
wave  on  a sharp  front , and  obtained  the  following  results: 

(1)  The  reflection  coefficient  for  energies  in  the  incident 
and  reflected  wave  pulses  is  er*/e^*  = (ti)j*/u>r*  )|  r|2  for 
the  oncoming  ionization  front,  and  = (u>r*/u)j,*  )|  rj 2 

for  the  recombination  front  case,  where  lr|2  = 1 when  the 
plasma  is  overdense  [cf.  Eqs.  (28)  and  (53)]. 
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(2)  The  duration  of  the  reflected  pulse  is  shortened  (lengthened) 
for  the  ionization  (recombination)  front  case,  so  that 


Tr*/Tj*  = uK*/<ur*,  and  the  power  reflection  coefficients  are 

P ,*/P  * = | Tj  2 for  the  ionization  front,  and  P */P.#  = 
r i r i 

(tu  */<i>i*)|  r|  2 for  the  recombination  front  [cf.  Eqs.  (29)  and 
(34)]. 

(3)  Energy  conservation  was  verified.  In  particular,  for  the 
ionization  front  in  the  overdense  case  ( | T | 2 =»  1 ) , it  was 
shown  that  the  energy  loss  upon  reflection  was  due  to  the 
excitation  of  a "magnetic  wave"  in  the  plasma  region.  The 
existence  of  this  wave  seems  to  have  gone  unrecognized 
until  now. 

We  next  considered  a front  with  an  arbitrary  electron  density 
profile.  We  found  that  conclusions  (l)  and  (2)  remain  valid  even 
in  this  more  general  situation.  Point  (3)  however  was  found  to  be 
altered: 

(J+)  As  the  thickness  of  the  front  increases  the  excitation 
coefficient  of  the  magnetic  wave  decreases  due  to  phase 
mixing  of  the  transverse  currents.  This  phase  mixing 
leads  to  a type  of  collisionless  heating  in  which  the 
transverse  temperature  of  the  electron  fluid  is  increased 
downstream  from  the  ionization  front. 

For  oblique  incidence  and  incidence  in  a waveguide  with  E transverse 
to  the  electron  density  gradient,  conclusion  (1)  was  again  found  to  hold. 
However,  for  the  waveguide  case  (2)  no  longer  applies,  since 
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Fig.  2 — Integration  path  C for  integrals  in  Eqs.  (58)  and  (59).  Poles  of  the  integrand 
of  (58)  are  denoted  by  x.  The  integrand  of  (59)  has  the  same  poles  as  that  of  (58)  but 
with  an  additional  pole  at  s * 1— koa  — iajc/U,  denoted  by  a dot  on  the  figure. 
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